Thermal phase transitions for Dicke-type models in the ultra-strong coupling limit 
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We consider the Dicke model in the ultra-strong coupling limit to investigate thermal phase 
transitions and their precursors at finite particle numbers A'' for bosonic and fermionic systems. We 
derive partition functions with degeneracy factors that account for the number of configurations 
and derive explicit expressions for the Landau free energy. This allows us to discuss the difference 
between the original Dicke (fermionic) and the bosonic case. We find a crossover between these two 
cases that shows up, e.g., in the specific heat. 
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The Dicke model, in its original weak-coupling and 
multi-mode form, l1| , has a long history as a paradigm for 
collective dissipation 01 , and, in its single mode form, as a 
test-bed for fundamental concepts such as the quantum- 
classical relation scaling [sj, or entanglement [oj 
near quantum phase transitions. Interest in the Dicke 
superradiance model has been furthered by the recent 
discovery of the 'Hepp-Lieb'-type quantum phase tran- 
sition p, 0| with Bose-Einstein condensates in an op- 
tical cavity Q. As a mean-field-type phase transition, 
the full phase diagram in the temperature-coupling con- 
stant plane was derived early pi]| - [l2| in the thermody- 
namic limit of — oo particles, cf. also [l^. In this 
Brief Report, we re-examine the thermal properties of (a 
somewhat generalized version of) this model with par- 
ticular emphasis on the influence of the quantum statis- 
tics on the quantum phase transition. We only consider 
the ultra-strong coupling limit between the atoms and 
the light (corresponding to the superradiant phase at low 
temperatures), but we distinguish between various cases 
of N bosons or N fermions distributed among Ng two- 
level sites. In particular, we derive simple expressions 
for the thermodynamic partition sums that can be used 
to easily calculate thermodynamic quantities such as the 
specific heat at finite particle number N, and to follow 
an interesting crossover between the case of N bosons on 
A^s = 1 site and the original Dicke case of N fermions on 
Ns = N sites . 

As a starting point we use the Dicke Hamiltonian with 
the single bosonic mode a, of frequency w. The angular 
momentum operators Jz, describe an ensemble of A^ 
two-level atoms with a level splitting ujq. The single mode 
Dicke Hamiltonian is 

n^OJa^a + ujQj^ + ^(a + a^)(J+ + J-). (1) 

Q^^Jy and (7 = 

(a^ — a) rotates and polaron-transforms the Hamil- 
tonian into T-L' = UHW with 



A unitary transformation with U = e"^"^ 
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The Hamiltonian "H' can be used as a starting point for 
a perturbation theory in wq, i.e. around the limit of very 
large coupling g ^ oo between the angular momentum 
and the photon mode. In this limit, the physics is then 
determined by a (trivial) free photon Hamiltonian ua^a 
and the angular momentum part oc J^. The analysis of 
the properties of the thermodynamic partition sum 
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for different physical realizations of Jz is the aim of this 
Brief Report. 

The role of particle statistics in can be qualita- 
tively understood by considering the existence or other- 
wise of a thermal phase transition for in the limit of 
A — )• oo. To this end, let us first recall that in the original 
Dicke model, Jz = ^ J2n=i '^n '^"^ the sum of A^ individ- 
ual (pseudo)-spin-i operators. Superradiant states with 
maximal spin polarization are then energetically favored 
by "Hjv, but there are only two configurations (all spins 
pointing either up or down) where that is achieved. All 
other spin configurations have larger (non-zero) entropy 
such that thermal fluctuations trigger a phase transition 
to a thermally disordered (normal) phase above a critical 
temperature. 

On the other hand, if the state space is restricted to the 
highly symmetric Dicke states | J, M) , the configuration 
space is much smaller and there is no gain in entropy for 
states with higher energy. In this situation, which corre- 
sponds to A^ bosons occupying either the upper or lower 
level of a two-level system, a thermal phase transition 
does not occur. 

The partition function and the number of configura- 
tions . - We start with discussing the degeneracy factors 
that appear when evaluating the partition sum Z^- Let 
us assume a configuration space with A^ particles dis- 
tributed among A^ two-level 'sites' i, all of which have 
the same up-level f and down-level energies. We write 
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such that 
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where c„ is the number of configurations with n particles 
in the down levels \, and — n in the up levels f- 

Fermions .- For fermions, the number of configura- 
tions for n particles in the down levels is (''^°), and for 
the remaiming N — n particles in the up levels it is 
{n-t)'^ consequently c„ = (^')(^'„)- n) For fermions 
with no two particles on the same site i, one has a re- 
stricted choice once all the f (or the \.) are occupied 

and thus c„ = Cn){^In)- ^'^^ ^^"^ particular case 
Ns = N, this corresponds to the N localized and dis- 
tinguishable (pseudo) spins in the original Dicke model 
discussed above and gives the partition sum 
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Bosons .- For N bosons in m single particle lev- 
els, there are (^^™r^) configurations which in our case 
means 



n + N,-l\/'N-n + Ns-l 



(6) 



The particular case of a single site Ng — 1 yields c„ = 1 
and thus 



Zbos 
N — 



N 

Tl = 
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(7) 



We also obtain this result by representing — ^{2b^b~ 
N) via one Holstein-PrimakofF boson or alternatively 
by writing Jz = ^(6|6-t- — ^Ifoj.) with two Schwinger bo- 
son modes for f and I by using -I- 6^6^ = N and 
noticing that the number n of J,-bosons uniquely fixes a 
configuration which means c„ = 1. 

Landau free energy function .- The most transpar- 
ent way to discuss the difference between the two cases 
Eq. ([5]) and Eq. ^ and generalizations thereof is by 
transforming the partition sums into integrals over an 
order parameter y. As we are dealing with a zero- 
dimensional field theory here, this is particularly sim- 
ple and is formally achieved by a Hubbard-Stratonovich 
transformation which here is simply given by the Gaus- 
sian integral identity, e^ = ^ dj/e^^'+^^'', applied 
to the respective Boltzmann factors in Zjv- For the dis- 
tinguishable (Dicke) case, we use the binomial formula 
as T,n=o (^)e^(^-2n) ^ (2coshz)^ to obtain (after sub- 
stituting y VN y) 



$°(2/) = y2_in(2cosh(2ay)), a = ff\/-, (8) 
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FIG. 1. The specific heat per particle (fcs = 1) as a function 
of Q^"^ = Ljj/{g'^l3) (proportional to the temperature) in the 
Dicke case of A*' localised spin-is, (derived from the fermionic 
partition sum Eq. (O), for iV = 20, 75, 400 and 1000 particles 
(dotted line; TV — >■ oo). Upper inset: shift of peak from 
critical point a* — l/y/2 for all A*' < 1000. Lower inset: peak 
height at maximum. Dashed lines: fit of the numerical data 
with /^pp,, = 1.64 - ^ and = 1.98 - 



where we introduced the dimensionless coupling param- 
eter a. 

Similarly, for the Ng = \ boson case, we carry out the 



geometric progression X]n=o 
l)/(e2^ - 1) to obtain 
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1. 



The so defined Landau free energy functions $(?/) now 
allow us to elucidate the critical properties of the models. 

First, we observe that in the (Dicke) case of distin- 
guishable particles, ^^{y) in Eq. ([5]) is A^- independent, 
and ^^{y) = (if is determined by the — 0~limit of 
the usual mean-field {N — >■ oo) expression for the free 
energy / per particle [Tol - [l^ . As expected, the term 
In (2 cosh (2a?/)) therefore is the mean field free energy 
of a single (pseudo) spin in the fluctuating field y. The 
parameter a determines the shape of ^^{y) and, in the 

— > oo limit, the position ?/d of the minimum of ^^{y) 
relevant for the asymptotic expression of the integral 
according to the Laplace method, as given by the self- 
consistent equation 



?/D — atanh2Q;?/D- 



(10) 



This has the unique solution j/d = when 2a^ < 1 cor- 
responding to temperatures T > Tc = 2g^ juj larger than 
the critical temperature Tf.. This solution describes the 
normal phase whereas for T < Tc there exist two minima 
in $^(2/) that describe the symmetry-broken superradi- 
ant phase. 

This is in contrast to the bosonic case, where the Lan- 
dau free energy 

^'^°'™(y) in Eq. ^ is iV-dependent 
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FIG. 2. The specific heat per particle as in Fig. HI) but for 
the bosonic partition sum Eq. ((7|), for — 2,5,12,20 and 
40 particles. Dashed lines: linear fit f^pp^r = 0.6 + 0.2A^ for 
maximum position (upper inset), fuiwer = 1^ fit for maxi- 
mum (lower inset). 



but acquires a simple form in the thermodynamic limit 
^ oo 

*!^°'°"(2/)-y'"2a|2/|, (11) 

which has two unique minima at y = ±a regardless of the 
value of a. This means that in the bosonic case, one is 
always in the superradiant, symmetry-broken phase and 
no thermal phase transition into a normal phase occurs. 

Thermodynamic behavior. - In the following, we discuss 
the specific heat C = P^d^ \og{Z]y) for both cases (we set 
fcs = 1). In the Dicke case (spin-^s corresponding to the 
fermionic partition sum Eq. ([5])), Fig. ([1]), the specific 
heat per particle can be calculated from 3>^(y) = /?/ in 
the thermodynamic limit A'^ — > cxd by eliminating deriva- 
tives oi yo and solving the self-consistent Eq. (fTOl) . lead- 
ing to the expected singular behaviour at the transition 
point a* = 1/V2. A 1/N expansion of Z^, Eq. (gl), 
works only well not too close to the critical point, as we 
have checked. Numerically, for finite N the peak height 
of C/N and the shift of the peak from a* give, however, 
a very good agreement with a 1/Vn correction fit. 

In the bosonic case. Fig. ([2]), the specific heat per par- 
ticle shows a totally different behavior. For small num- 
bers of bosons {N = 5, 6, ...12), C/N first increases with 
N and then (for N > 12) decreases with as does its 
maximum. The position of the maximum of C/N is 
linear the particle number N. In the limit N ^ oo the 
specific heat per particle is zero which can be understood 
from Eq. ()11|) by using the above mentioned thermody- 
namic relations and considering that y — ±a. 

Crossover between bosonic and Dicke case We now 
discuss an interesting crossover between the two cases ob- 
tained above by regarding the bosonic case with A^^ > 1. 
For Ng of the order of the particle number TV, we expect 
the bosons to spread over many energetically equivalent 
configurations which, due to the entropy argument given 
in the introduction, should re-establish the thermal phase 
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FIG. 3. Crossover between bosonic and Dicke case in the 
Landau free energy, Eq. (|12|) for A'^ = 400 bosons and coupling 
strengths a — 0.2, 0.4, 0.6, 1.0 within each graph. Upper left 
Ns = 1 {Na'. number of two- level sites), this corresponds to 
the simple bosonic case Eq. ((9]) (and Eq. (|11|) for — >■ oo), 
i.e. the absence of a phase transition, since the number of 
minima remains two irrespective of a. The phase transition is 
re-established for macroscopic degeneracy (upper right Ns = 
A^/4, lower left Ns — N), since a transition from a double- 
to single-minimum occurs. The lower right shows the Dicke 
case ^^{y), Eq. ([H]), with the phase transition occurring at 
a* = 1/V^Ri 0.71. 

transition found in the Dicke (fermionic spin-^) case. 
We therefore generalize the Landau free energy for the 
bosonic case Eq. ([9]) to arbitrary Ng, 

1 ^ 

$^X(y) ^y'-j^ In^ c„e2"''(^-2") (12) 

with c„ given by Eq. ([6]). Fig. (|3]) shows the crossover 
in the free energy when passing from the bosonic case 
with Ns — 1 (showing no phase transition when vary- 
ing a) to larger degeneracies, where for Ng ^ A'', the 
free energy <&^'^"(y) becomes equivalent to the one of 
the Dicke case, Eq. This can be understood by using 
Stirling's formula to expand the number of configurations 
c„, Eq. dSJ, which leads to c„ - {^){Ns/N)^ and there- 
fore the partition sum 

- (j^) Z^, Ns » A^ » 1. (13) 

Correspondingly, in this limit the bosonic free en- 
ergy per particle differs from the Dicke free energy by 
— kBT\n{Ns/N) which just describes an additional en- 
tropy term caused by the enhanced 'volume' of configu- 
rations. Thus for Ns ^ A^, multiply occupied sites play 
no role any longer, nor does the statistics of the particles 
(fermions or bosons), and the only remaining relevant 
statistical information is the number of occupied upper 
and lower levels as in the Dicke case. As the entropy 
gain kB^n{Ns/N) is just a constant, the specific heat 
C/N then has to coincide with the specific heat in the 
Dicke case. 
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FIG. 4. The specific heat per particle at fixed particle num- 
ber = 400 for different number of two-level sites Ns — 1 
(bosonic case, bottom line), A^s = 100, Ns = 400 (Dicke 
case) and Ns = 600,5000 as obtained from Eq. ©. Up- 
per inset: with increasing A'^a the peak position in C/N 
shifts towards the critical point = 2 of the Dicke case. 
Lower inset: increase of maximum of C/N with Ng. Dashed 
lines: fit of the numerical data with /upper = 1-82 — ^ and 

r4 — l.SiVa 

Jlower 458.6-l-iVs ' 

This behavior of the specific heat C/N is shown in 
Fig. (jH), where we vary Ng from the bosonic case A^s = 1 
over Ng = N to finally higher values of iV^ ^ TV cor- 
responding to the Dicke (fermionic) form of C/N from 
Fig. (P). With N fixed, the C/N curves have peaks that 
shift from larger to smaller values of a"^ with increasing 
Ng/N corresponding to a decreasing N/Ng in the A^^ = 1 
bosonic case in Fig. At even larger Ng/N ^ 1 the 
curves approach the Dicke form at the critical point a* of 



the Dicke thermal phase transition. Note that the peaks 
approach a* from the right here, whereas in the Dicke 
case. Fig. ([1]), Ng = N was fixed and the C/N peaks 
approached a* from smaller with increasing N. 

Conclusion and outlook . - Our results indicate that the 
Dicke model displays an interesting thermodynamic be- 
havior if one considers the possible configurations of N 
particles in A^^ two-level systems. In the ultra-strong 
coupling limit, the model reduces to an atomic self- 
interaction term for which we have derived explicit ex- 
pressions for the partition function. Landau free energy 
and specific heat. We find a crossover in the specific 
heat from a bosonic form at small Ng (where no thermal 
phase transition occurs) to the form following from the 
(original) Dicke model for Af^ ^ iV ^ 1 that displays a 
thermal phase transition at a^^ = kBTuj/g^ = 2 between 
the normal and the superradiant phase. 

In our calculations, we regarded Hn in Eq. ^ as an ef- 
fective Hamiltonian for the ultra-strong coupling regime 
between atoms and light of the Dicke model Eq. (jlj. 
In analogy with the simulation of the Dicke model with 
Bose-Einstein condensates Q , we suggest the various lim- 
its of Hn to be regarded as effective models to be sim- 
ulated with, e.g., cold atoms. A further challenge would 
be an extension to Dicke- type models with finite coupling 
strengths, multiple levels [IJI, and degeneracies Ng in the 
bosonic case. 
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